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ABSTRACT 

We link the evolution of the galaxies in the hierarchical clustering scenario with the changing accretion 
rates of cold gas onto the central massive black holes that power the quasars. We base on galaxy 
interactions as main triggers of accretion; the related scaling laws are taken up from Cavaliere & Vittorini 
(2000), and grafted to a semi-analytic code for galaxy formation. As a result, at high z the protogalaxies 
grow rapidly by hierarchical merging; meanwhile, much fresh gas is imported and also destabilized, so 
the holes are fueled at their full Eddington rates. At lower z the galactic dynamical events are mostly 
encounters in hierarchically growing groups; now the refueling peters out, as the residual gas is exhausted 
while the destabilizing encounters dwindle. So, with no parameter tuning other than needed for stellar 
observables, our model uniquely produces at z > 3 a rise, and at z < 2.5 a decline of the bright quasar 
population as steep as observed. In addition, our results closely fit the observed luminosity functions of 
quasars, their space density at different magnitudes from z « 5 to z w 0, and the local ttibh — & relation. 

Subject headings: galaxies: active - galaxies: formation - galaxies: evolution - galaxies: interactions - 
quasars: general 



1. INTRODUCTION 



Observations pinpoint massive dark objects, conceiv- 
■£3d>ly supermassive black holes (BHs) with masses itlbh ~ 
^0'' — 5 10 9 Mq, at the center of most nearby bright galax- 
• les (see Richstone et al. 1998). On the other hand, the 
^space density of bright optical quasars (QSOs) is found to 
•JJpe lower by some 10 -2 relative to the bright galaxies (see, 
^.g., Boyle et al. 2000). 

03 These observations support the view of the QSOs as a 
short active phase (lasting At ~ 10 8 yrs) of supermassive 
BHs which accrete surrounding gas at rates rh acc ~ 1 — 10 2 
M Q /yr, see Rees (1984). If the accretion history starts 
from small primordial seeds (Madau & Rees 2000), the ac- 
cretion rate sets not only the QSO bolometric luminosities 
L = rjc 2 m acc <; 10 48 erg/s (with standard mass-to-energy 
conversion efficiency r/ w 0.1, see, e.g., Yu & Tremaine 
2002), but also the relic BH masses. 

Additional observations relate the history of m acc (z), 
and so the QSO evolution, with the galaxy structure and 
formation. These include the very magnitude of m acc that 
needs to be drawn from the whole bulge of the hosts, and 
the correlation of tubh with the bulge luminosities (see 
Richstone et al. 1998), or better with the velocity disper- 
sion a (Ferrarese & Merritt 2000; Gebhardt et al. 2000). 

Relations of the sort have been explored in the frame- 
work of the hierarchical clustering scenario; here the galax- 
ies condense from gravitationally unstable regions in the 
primordial density field of the dark matter (DM), then 
are assembled into poor groups and eventually into rich 



clusters. Such attempts are based either on analytic ap- 
proaches (see Cavaliere & Szalay 1986; Efstathiou & Rees 
1988; Carlberg 1990; Haehnelt & Rees 1993; Haiman & 
Loeb 1998; Burkert & Silk 2001; Hatzminaoglou et al. 
2002), or on semi-analytic models (SAMs) of galaxy forma- 
tion (see Kauffmann & Haehnelt 2000; Monaco, Salucci, 
Danese 2000; Volonteri, Haardt & Madau 2002). But the 
fraction / = Am acc /m c actually accreted out of the galac- 
tic gas remained as a phenomenological quantity, described 
only in terms of scaling recipes, either containing tun- 
able parameters or adjusted from the outset to the form 
of the observed correlations. Yet, all such models hardly 
accounted for the observed dramatic drop of the QSO pop- 
ulation between z « 2.5 and the present. 

On the other hand, a physical law for the accretion has 
been proposed by Cavaliere & Vittorini (2000, CV00); they 
suggested the amount of the cool galactic baryons accreted 
onto BHs to be given by the fraction / destabilized and 
sent toward the nucleus by the gravitational torques aris- 
ing cither in major merging events (at z J> 3) or in en- 
counters of the host with other galaxies inside a common 
DM halo (at z <J 2.5). They use two analytic evaluations 
for the two regimes; in the latter they derive declining 
encounter rates and fast gas exhaustion concurring to pro- 
duce an interestingly steep evolution of bright QSOs. 

Here we self-consistently compute halo and galactic 
quantities using a hierachical SAM of the kind proposed 
by Kauffmann et al. (1993); Somerville & Primack (1999); 
Cole et al. (2000). Specifically, we insert the CV00 model 
in the SAM developed by Menci et al. (2002); this in- 
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eludes galaxy interactions, and accounts for several sta- 
tistical properties of the galaxy population like counts, z- 
resolved luminosity functions (LFs), sizes. This will allow 
us to continuously link the accretion history rh acc (z) and 
the QSO evolution with a fiducial model of galaxy forma- 
tion consistent with the stellar observables. 

2. MODELING THE GALAXY EVOLUTION 

The independent variables of our SAM (Menci et al. 
2002) arc: the circular velocity V of the host DM halos 
(groups and clusters of galaxies with mass M and virial 
radius R) containing the galaxies; the circular velocity v 
of the DM clumps associated with the individual mem- 
ber galaxies. The former grow hierarchically to larger 
sizes through repeated merging events (at the rate given 
in Lacey & Cole 1996), while the latter may coalesce ei- 
ther with the central galaxy in the common halo due to 
dynamical friction, or with other satellite galaxies through 
binary aggregations. The timcscalc for such galactic pro- 
cesses generally exceed the timescale for the merging of 
the host halos, so member galaxies accumulate in growing 
host halos. 

At the cosmic time t the SAM yields the following dy- 
namical quantities. 

• The distribution N(v,V,t) (per Mpc 3 ) of galaxies with 
circular velocity in the range v — v + dv in host halos with 
circular velocity in the range V — V + dV. This is de- 
rived by computing iteratively the probability for a host 
halo to be formed from its progenitors, together with the 
probability that the member galaxies coalesce due to ei- 
ther dynamical friction or binary aggregations. As initial 
condition we assume a Press & Schechter distribution, and 
we assign one galaxy to each host halo. From N(v,V,t) 
we derive the number NT(V,t) of galaxies in a host halo 
(membership), and the overall distribution of galaxy cir- 
cular velocity N(v,t) irrespective of the host halo. 

• The tidal radius r t {v), and the disk radius rd{v) and 
rotation velocity Vd(v) computed after Mo, Mao & White 
(1998). 

• The average relative velocity V re i(V) of the galaxies in 
a common DM halo. 

These dynamical quantities are associated to the 
baryons contained in the galactic DM halos in the way 
standard for SAMs. Initially the baryons contained in a 
galactic halo are in the amount mflb/Q (with man 3 
the DM mass of the galaxies) and at the virial tempera- 
ture; the mass m c of cold baryons is that residing in re- 
gions interior to the "cooling radius" . Stars are allowed to 
form with rate to* = (m c /td yn ) (f/200kms -1 ) - "*, with 
the disk dynamical time evaluated as td = rd/vd- Finally, 
a mass Ato/j = (3 to* is returned from the cool to the hot 
gas phase due to the energy fed back by canonical type 
II Supernovae associated to to*. The feedback efficiency 
is taken to be (3 — (v/vh) ah - The values adopted for the 
parameters a* = —1.5, ah = 2 and Vh = 150 km/s fit both 
the local B-band galaxy luminosity function and the Tully- 
Fishcr relation, as illustrated by Menci et al. (2002). At 
each merging event, the masses of the different baryonic 
phases are refueled by those in the merging partner; the 
further increments Am c , Am*, Am/, from cooling, star 
formation and feedback are recomputed iterating the pro- 
cedure described above. 

All computations are made in a A-CDM cosmology with 



Q,q = 0.3, Q\ — 0.7, a baryon fraction fi& = 0.03, and Hub- 
ble constant h = 0.7 in units of 100 km s _1 Mpc -1 . 

3. BH ACCRETION TRIGGERED BY GALAXY 
ENCOUNTERS 

Here we just recall the key points of CV00, and recast 
their main equations in a form suitable for implementation 
in our SAM. 

For a galaxy with given v, the accretion of a fraction 
/ of the cold gas onto the central BH is intermittently 
triggered by interactions occurring at a rate 



n T {V)H{v,V)V rel (V) 



(1) 



Here tit — Nt / (iHR 3 / 3) , and the cross section T,(v, V) w 
7r( (rf + r t 2 )) is averaged over all partners with tidal ra- 
dius r' t in the same halo V. To this effect, the sizes, the 
relative velocities of the galaxies and the distribution of 
galaxy circular velocities in a common halo are previously 
computed from the SAM sector described in §2. 

The fraction of cold gas accreted by the BH in each in- 
teraction event is computed in eq. A3 of CV00 in terms the 
variation Aj of the specific angular momentum j w Gm/vd 
of the gas, to read 



f(v,V) 



Aj 
3 
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to b VI 



(2) 



Here b = max[rd, d(V)\ is the impact parameter, evaluated 
as the maximum between the galaxy and the average 

1/3 

distance d(V) = R/N T (V) of the galaxies in the halo. 
Also, to' is the mass of the partner galaxy in the interac- 
tion, and the average runs over the probability of finding a 
galaxy with mass to' in the same halo V where the galaxy 
to is located. The prefactor accounts for the probabil- 
ity (1/2) of inflow rather than outflow related to the sign 
of Aj, and for the fraction (1/4) of the inflow actually 
reaching the nucleus rather than kindling circumnuclear 
starbursts, see Mirabel & Sanders (1996). 

Note from cqs. (1) and (2) that both the interaction 
rate r" 1 and the accreted fraction / decrease with time, 
since in the growing host halos the increasing membership 
Nt(V) is offset by the increasing R, V and V re i. In a 
group the values for / range from a few 10 _1 to several 
10~ 3 . The above behaviour of the accreted fraction / is 
illustrated in the top panel of fig. 1. 

The average gas accretion rate triggered by interactions 
at z is given by (CV00 eq. 5) 



m a 



z(v,z) 



/ f{v,V) m c {v) 

\ T r (v,V) 



(3) 



where the average is over all host halos with circular ve- 
locity V. The bolomctric luminosity so produced by the 
QSO hosted in a given galaxy is then given by 



L(v,t) 



r]c 2 Am 

a 



(4) 



Here t w td ~ 5 10 7 it/to) yrs is the duration of the ac- 
cretion episode, i.e., the timescale for the QSO to shine; 
ATO acc is the gas accreted at the rate given by eq. (3). We 
take r] = 0.1, and we obtain the blue luminosity Lb by 
applying a bolomctric correction of 13 (Elvis et al. 1994). 
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The mass of the BH hosted in a galaxy with given v at 
time t is updated after 

rriBH{v,t) = (1 - 77) / m acc {v,t')dt' (5) 
Jo 

assuming in all galaxies small seeds BHs of mass 10 2 M Q 
(Madau & Rees 2001); our results are insensitive to the 
specific value as long as it is smaller than 10 5 M Q . 
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Fig. 1. - Top: The fraction of cold gas accreted onto BHs as a func- 
tion of the circular velocity of the host galaxy at z = 0.4 (solid line), 
z = 1 (dotted line) and 2 = 3 (dashed line). The inset shows the 
the cosmological density of all the available cold gas in galaxies as 
a function of the redshift z. Bottom: The mass function of BHs for 
various values of the redshift z. 

4. RESULTS 

In fig. 1 (bottom) we plot the mass function of BHs 
at four different redshifts. Note the mild evolution from 
z = 1 to the present, due to the decrease of the accreted 
fraction / (see top of fig. 1) and of the related Am acc dut 
to the declining rate of merging and encounters (see §3), 
and to the simultaneous exhaustion of the galactic cold 
gas m c {z) available for accretion (see inset in fig.l). 

The relation we obtain between the BH mass and the 
central 1-D velocity dispersion a of the host galaxy is 
shown in fig. 2. We use the canonical relation a = v/y/2 
approximately holding for an isothermal profile (see Bin- 
ney & Tremaine 1986). 

The steep slope of our relation (close to Mbh oc a 4 in 
the range 10 7 M Q <; Mbh <, 10 9 M ) is the combined re- 
sult of two processes: the merging histories of the galactic 
DM clumps, which by themselves would imply the mass 



of cold available gas to scale as a 2 5 ; the destabilization of 
the cold gas by the interactions, which steepens the rela- 
tion to a 3 5 . The further steepening to cr 4 is provided by 
the Supernovae feedback which depletes the residual gas 
content in shallow potential wells. No effort has been done 
to adjust the latter to match the relation to the data; a 
slope steeper yet would be provided, especially in the up- 
per range, by adding the feedback from the QSO emission 
itself onto the host gas (which does not alter the average 
galactic star formation in a group). 
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Fig. 2. The derived relation between the BH mass and the 1-D 
velocity dispersion of the host galaxy (solid line) is compared with 
data from Ferrarese & Merritt (2000, filled squares) and Gebhardt 
et al. (2000, circles). 

The evolving LF of the QSOs is derived from N(v, t) by 
applying the appropriate Jacobian, see §3. The LF will 
include a factor r/(r r (u)) < 1 since the luminosities in eq. 
(5) last for a time r and are rekindled after an average 
time r r . The result is 



N(L,t)=N(v,t) T - 

(T r ) 



dv 
dL 



(6) 



and is shown in fig. 3. Note the strong evolution from 
z k 0.5 to z k 2, due to both the declining Am acc (z) and 
the timescale r shortening at higher z. 

In fig. 4 we plot the predicted cosmic density pq(z) of 
QSOs brighter than M B = -24 and M B = -26. The com- 
puted emissions, limited by the Eddington luminosity Le, 
peak at z i=s 2.5; such redshifts are those where our SAM 
predicts the buildup of bright galaxy to peak (see Menci 
et al. 2002). Note that if no Eddinton limit were imposed, 
the density would be that represented by the dashed line 
in the bottom panel, which overshoots the data; this shows 
that at z > 2.5 the QSOs do emit at their full Le, due to 
the high merging rate which provides abundant refueling. 
At later times, the fueling is triggered mostly by the inter- 
actions with encounter rate declining as discussed in §3, 



4 



QSO Evolution in Hierarchical Structures 



and draws from scarce residual gas approaching exhaus- 
tion. This causes the dramatic fall of the QSO density we 
predict for redshifts z < 2. 
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Fig. 3. The LFs from our model (solid lines) are shown for z = 0.55 
(lower curve), z = 1.2 (middle curve) and z = 2.2 (uppermost 
curve), and are compared with the data points. These arc taken 
from Hartwick & Shade (1990, solid squares) and Boyle et al. (2000, 
crosses), and rescaled to our cosmology with Qo = 0.3, Q\ = 0.7, 
h = 0.7. We also show as a dashed line the model LF for z = 4.2 
compared to the Sloan data from Fan et al. (2001, diamonds). The 
dotted line is the predicted LF for z = 3.4. 

5. CONCLUSIONS 

Our approach links the evolution of QSOs to that of 
their host galaxies on using scaling laws for the BH accre- 
tion based on triggering galaxy interactions, with no new 
parameters added to the galaxy formation model. It pro- 
vides good fits to the observations concerning the itibh — v 
correlation, the evolution of the QSO luminosity function, 
and their space density selected at different optical magni- 
tudes. The agreement holds in the full range z < 6 contin- 
uously covered by our model, and supports the underlying 
picture that follows. 

At z > 2.5 the QSOs emit at their full Eddington limit, 
because the high merging rate in this epoch of galaxy as- 
semblage insures both a high baryonic content in their 
hosts and an abundant BH fueling. As shown in fig. 3, in 
this range the LFs are flat (consistent with by recent data, 
see Fan et al. 2003) and rise with time. At z < 2.5 the 
steep drop we find for the population results from three 



concurring processes: 1) the declining rate of merging be- 
tween early galaxies, which halts the acquisition of new 
gas available for accretion; 2) the progressive exhaustion 
of the baryon reservoirs in the hosts, consumed by fast 
conversion into stars and by previous accretion episodes; 
3) the eventual decline of the fraction / ~ Aj/j of residual 
cold gas which is destabilized and accreted onto the cen- 
tral BHs by the dwindling interactions between galaxies. 
The latter (see fig.l) is at variance with the Kauffmann 
& Haehnel (2000) model, and is crucial in matching the 
observed strong evolution of QSOs at z < 2; as a con- 
sequence, our model predicts an average Eddington ratio 
dropping from L/Le ~ 1 at z w 2.5 to L/Le ~ 10~ 2 at 
z w 0, with a weak dependence on rriBH, consistent the 
with available data (see Woo & Urry 2002). 
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Fig. 4. The predicted cosmic density of bright (upper) and very 
bright (bottom panel) QSOs. The dashed line represents the out- 
come when no Eddington limit were assumed (see text). The den- 
sities have been rescaled as appropriate for the critical cosmology 
with h = 0.5 used by the authors for their (magnitude-limited) data; 
these are taken from Hartwick & Shade (1990, solid squares); Gold- 
schmidt & Miller (1998, triangles); Warren, Hewett & Osmer (1994, 
open circles); Schmidt, Schneider & Gunn (1995, filled circles); Fan 
et al. (2001, crosses). 

The above picture implies a specific connection of the 
QSO emission with observable properties (magnitudes, 
colors and star content) of their host galaxies. On this 
we shall present elsewhere our results; we only stress that 
in our present model the host galaxies of bright QSOs have 
dynamic properties conducive to strong interactions, and 
hence to high luminosities, at early z. 



REFERENCES 



Boyle, B.J., Shanks, T., Croom, S.M., Smith, R.J., Miller, L., 

Loaring, N., Heymans, C. 2000, MNRAS, 317, 1014 
Burkert, A., & Silk, J., 2001, ApJ, 554, L151 
Carlberg, R.G., 1990, ApJ, 350, 505 
Cavaliere, A., & Szalay, A., 1986, ApJ, 311, 589 
Cavaliere, A., Vittorini, V., 2000, ApJ, 543, 599 



Cole, S., Lacey, C.G., Baugh, CM., Frenk, OS., 2000, MNRAS, 319, 
168 

Efstathiou, C, & Rees, M.J., 1988MNRAS, 230, 5 
Elvis M. et al. 1994, ApJS, 95, 68 

Fan, X. et al., 2001, ApJ, 121J54 

Fan, X. et al., 2003, preprint jastro-p h/0301135 



Menci et al. 



5 



Ferrarese, L. Merritt, D., 2000, ApJ, 539, L9 
Gebhardt, K. et al., 2000, ApJ, 539, L13 
Goldschmidt, P., Miller, L., 1998, MNRAS, 293, 107 
Haehnelt, M., & Rees, M.J., 1993, MNRAS, 263, 168 
Haiman. Z., & Loeb, A., 1998, ApJ, 503, 505 

Hatzminaoglou, E., Mathez, G., Solanes. J.-M.. Manrique, A., 

Salvador-Sole', E., 2002, prpeprint astro-ph/0212002 
Hartwick, F.D.A., Shade, D., 1990, AHa&a, zs, 43Y 
KaufTmann, G., White, S.D.M., & Guiderdoni, B., 1993, MNRAS, 
264, 201 

KaufTmann, G., Haehnelt, M., 2000, MNRAS, 311, 576 
Lacey, C, & Cole, S., 1993, MNRAS, 262, 627 
Madau, P., & Rees, M.J., 2000, ApJ, 551, L27 



Magorrian, J., et al. 1998, AJ, 115, 2285 

Menci, N., Cavaliere, A., Fontana, A., Giallongo, E., Poli, F., 2002, 
ApJ, 578, 18 

Mirabel, I.F., Sanders, D.B., 1996, ARA&A, 34, 749 

Mo, H.J., Mao, S., & White, S.D.M., 1998, MNRAS, 295, 319 

Monaco, P., Salucci, P., Danese, L., 2000, 311, 279 

Rees, M.J., 1984, ARA&A, 22, 471 

Richstone, D., et al. 1998, Nature, 395, 14 

Schmidt, M., Schneider, D.P., Gunn, J.E., 1995, AJ, 110, 68 

Somerville, R.S., & Primack, J.R., 1999, MNRAS, 310, 1087 

Yu, Q., & Tremaine, S., 2002, MNRAS, 335, 965 

Warren, S.J., Hewett, P.C., Osmer, P.S., 1994, ApJ, 421, 412 

Woo, J.-H., & Hurry, CM., 2002, ApJ, 579, 530 



